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Abstract 



We have discussed earlier the correlation functions of the random vari- 
ables det(A — X) in which X is a random matrix. In particular the moments 
of the distribution of these random variables are universal functions, when 
measured in the appropriate units of the level spacing. When the A's, instead 
of belonging to the bulk of the spectrum, approach the edge, a cross-over 
takes place to an Airy or to a Bessel problem, and we consider here these 
modified classes of universality. 

Furthermore, when an external matrix source is added to the probability 
distribution of X, various new phenomenons may occur and one can tune 
the spectrum of this source matrix to new critical points. Again there are 
remarkably simple formulae for arbitrary source matrices, which allow us 
to compute the moments of the characteristic polynomials in these cases as 
well. 
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1 Introduction 



In the theory of random matrices , the n-point correlation functions of the eigenvalues 
are known to be expressible as the determinant of a two-point kernel |jl|, |^. The 
expressions for those kernels depend on the various classes of universality : it is a 
simple sine-kernel within the bulk of unitary invariant ensembles, an Airy kernel at 
the edge of the spectrum or a Bessel kernel for other invariance properties of the 
measure. The level spacing probability p{s), has also been computed recently for 
those different kernels P, H- 

Another interesting object is given by the average moments of the characteris- 
tic polynomial of the random matrix. These characteristic polynomials have been 
first investigated in |^ for a uniform probability measure on unitary matrices, 
in connection with the moments of the Riemann zeta-function. These results have 
been generalized to random hermitian N x N matrices X with a unitary invariant 
probability measure 



Explicit formulae for the 2i^-point functions 

2K 

F2k{\i, ■ ■ ■ , X2k) =< n det(Ai -X)> (2) 

1 

have been derived, which show that these functions are universal in the Dyson limit, 
in which the size N of the matrices goes to infinity, the distances between the A's go 
to zero, and the products N{Xi — Xj) remain finite. In particular the moments 



of the distribution of the characteristic polynomials were given in the large N limit 




(1) 



F2k{X, ■ ■ ■ 



A) =< [det{\ -M)Y'' > 



(3) 



by S0 



exp -{NKV{X))F2k{X, ■ ■ ■ , A) = (27riVp(A))^'e-^^7^, 



(4) 



1 



with 

K-l ^, 

provided A belongs to the bulk of the support of the distribution of the eigenvalues, 
i.e. provided p(A) does not vanish. Then one sees explicitely that the only dependence 
upon the probability measure is through the average density of eigenvalues p(A), and 
even the coefficient is a universal number. 

However the result does take different forms for different universality classes. Our 
previous investigations for the three classical Lie groups, U(N), Sp(N) and 0(N), are 
extended here to the Bessel kernel and Airy kernel, for which the density of states p(A) 
presents a singularity at the edge of the spectrum. Furthermore we have considered 
a Gaussian case in which an external matrix source is present Q in the probability 
distribution of the matrix 

P{X) = - exp (-A^Tr-X^ + NTrAX). (6) 

Explicit and simple formulae will be derived here again for the correlation functions 
and the moments of the characteristic polynomials of the matrix X, which depend 
on the spectrum of the matrix A. By tuning the spectrum of A appropriately, one 
can generate a number of different situations. For instance, we have investigated in 
the past the case in which the average spectrum of X presents a gap in the presence 
of A, and by tuning A one can study the critical point at which this gap vanishes. 
This creates again a new class of universality, and a new kernel [^, TU|. Other cases. 



such as 2D gravity in the double scaling limit, or the Penner model, would certainly 
be of interest as well. 

2 Sine-kernel 

For completeness, and for later use, we begin with the bulk unitary case, governed 
by the sine- kernel, but with a derivation which differs from our previous one [^]. 
An interesting geometric interpretation of this problem will also be provided. The 



kernel, from which all the correlation functions may be obtained, is given in terms 
of orthogonal polynomials for finite N, but reduces in the Dyson large N-limit to the 
sine-kernel 

sin(3: - y) 

K{x,y) = 7 

x-y 

in which x and y are the eigenvalues measured in the scale of the average spacing 
(27rp(A)iV)"'. Then, one obtains the normalized moments 



^NKVix) <[det(A-X)p^> detir(A„A,-j 

where A(A) = ni<j(Ai — Xj) and i,j = 1, ...,K. The r.h.s. may be expressed as a 
contour integral, following eq.(52) of 0, 



This may be further reduced to 



Ik = det(a„m) (10) 



1 9" sm{u — v) 
ralm! aw"- af ™ u — v 
where n, m = 0, 1, i^T — 1. The explicit evaluation of the determinant of a^^rn gives 

det(a_) = 2^^-^ n'(;^^ (12) 

We do recover in this way the factor '^k (0)(up to a factor 2^^~^ due to a different 
normalization normalization of the kernel). 

It is quite remarkable that this universal normalizing factor 'Jk has a geometric in- 
terpretation as a Fredholm determinant of the Dirac Laplacian on the two dimensional 
sphere 5*^. The determinant of the Laplacian has been discussed in the connection 
to string theory |TT], |T2|, and the relation of 'Jk to this Fredholm determinant of the 



Laplacian has been noticed in 0. Indeed let us show that 



where 7 is Euler's constant and A^{z) the determinant of a Dirac operator, defined 
below. The derivation goes as follows : let us introduce a function G{z) which satisfies 
the functional relation 

G{z + l) = T{z)G{z). (14) 
It is then straightforward to verify that 



7k 



n 



K-l 



^K-2K 



2 7r^+^ 



G{K + 



(15) 



A function G, satisfying the functional relation (0), is known in the literature as a 
Barnes function (or as the inverse of a di-gamma function). It is defined by 
1 



G{z + l) 



(27r)^/^ e"M^+(i+'>')^'] n (1 + i)"e~'^+'='/(2") 

1 L ^ 



(16) 



It has been noticed earlier (|T3[)that this Barnes function is related to the Fred- 
holm determinant of the Laplacian on S'^ . Indeed this Fredholm determinant is the 
(regularized) product 

A(^) = n(i - (17) 

where the A/ are the eigenvalues of the Laplacian, and gi their degeneracy, i.e. A; = 
/(/ + 1) with multiplicity (7^ = 2/ + 1, / = 0, 1, 2, .... It is convenient to shift 2; by 1/4, 
since this yields the the spectrum of the Dirac operator 



Then the regularized (shifted) Fredholm determinant 



'1^ 



{I 



)e(lTT7iF] 



7l2i+l 



factorizes as 



A( V) = A+(zy)A+(-zy) 
with the determinant of the Dirac operator A+(z) given by [p 



A+(^) 



n[(i 

i=0 



_^gI+f72 + 2(i/l/2)^ j2i+l_ 



(19) 



(20) 



(21) 



Then this Dirac determinant A"*" is related to the Barnes function by 



Ta-z)Ga 



We thereby recover the expression relating the moment 7;^ to the determinant (p!3D. 

This relation between the moments of the distribution and the determinant of the 
Dirac operator on S"^ is in fact general. For instance in the simplest case of a single 
Gaussian random variable, the moments are 



|2 



oo ^ ]^ 

oo 2 



x2^e-"'dx = r(ir+^); (23) 



T{K + 1/2) is thus the equivalent of for this trivial problem. If we consider the 
"Laplacian", i.e. the harmonic oscillator whose eigenvalues are A„ = n, then the 
Fredholm determinant A (A) is 



oo \ 

A(A) = -An(l--)e^ 

n=l 



n 



Hence, we have 



r(-A) 



.2K 



(24) 



< > A(A) l^=-(^+^) ^^^^ 

The expression (^) is a multi-variable version of this Gaussian integral. 

An additional point of interest is that the Fredholm determinant of this Laplacian 
on 5*^ may be factorized further into a product of two factors ; it turns out that each 
factor enters into the corresponding expression for the symplectic and orthogonal 
cases respectively. This will be seen below when we examine the moments related to 
the Bessel kernel. 

3 Bessel kernel 



We have discussed in our previous work the ensembles invariant under the unitary 
symplectic and unitary orthogonal Lie groups 0. The kernels for those ensembles 

5 



are |T|, |T5|, |T6 

K(,,y) = LCjl^t:ZJ!l^^^^^±Jil) (26) 

ZTT X — y X + y 

where the minus sign corresponds to the Sp and the plus sign to the O ensemble. It 
is convenient to introduce the Bessel kernel defined by 

K.(x,y)=^"^^^^"^^^~^"^^^'^-^^^ (27) 

x-y 



Since Ji/2ix) = y2/7rxsinx, J-i/2ix) = y2/7rxcosa;, the kernels for the Sp and O 
ensembles are both related to this Bessel kernel 

K±{x,y) = ^K^y2{x\y^) (28) 

namely, a = 1/2 and a = —1/2 represent respectively the Sp and the O ensemble. 
We consider from now on an arbitary a. The 2K-th moment at the origin (A = 0) is 
expressed as 

= / / {2nyi]l,!^^vf^ n(^^^^0"^«(«^' (29) 
We define now the two functions (j){z) and ip{z) by 

UV~z) = i^TT^r^M (30) 

2 T{a + 1) 

and 

Their expansions in powers of x are given by 

Keeping aside trivial factors we are then led to the kernel Ka{x, y) defined as 

Kc^ix^y) = 2^^_y^ \^i^)i'iy) - i'{^)^{y)] (34) 
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As before, we have 



with 



Ik = det(a, 
1 9" d 



a. 



KJu^v) 



\u=v=0 



IK 



-K^-aK 



2K-1 

n 



(35) 
(36) 

(37) 



n\m\ du"- dv"^ 

This determinant may be computed exphcitly, and it is given by 

1 

1=0 (« + 0! 

(We have Ji = for a = 0, | and — |, respectively.) 

It is interesting to relate the three determinants that we have introduced hereabove 
for the sine-kernel and for the Sp and O cases. The determinant for the sine-kernel 

(in is 



det 



/ 1 






1 

3 




V 



In the symplectic case, a 



|, we have 



l-Sp 



det 



V 



'30 




J_ 

20 





20 
7! 





_ J_ 

30 



20 
7! 



...\ 
■■■) 



(38) 



(39) 



In the orthogonal case, the determinant becomes for a = — | 



2' 



/ 1 



20 



Iq = det 

V... 

Thus, we find the factorization of ( |55| ) as the product of 
numerical factor due to the normalizations, 



Isp X I 



o 



(40) 

and (^) , up to a trivial 
(41) 



The factors 7^- for the unitary, symplectic and orthogonal case are related 

7S?^^r, and = (HET^ /!)V(nlV^ ^!), 7?^^ = 2^^^+^)/^ D^i /!/ nf=i(2/)! and 
^{o)2i^(i^+i)/2-i-Qi^-i/i/-Q^^i(^2/)!. It is again remarkable that, for arbitrary a, 'Jk 
may still be expressed as the Fredholm determinant of the Laplacian, in which the 
eigenvalues are shifted by the amount a |T^. 



4 Airy kernel 



When the eigenvalues he near an edge Ac of the support of the asymptotic density of 
states (an edge of Wigner's semi-circle in the Gaussian case), in a neighbourhood of 
size iV"^/^ of that edge, there is a cross-over from the sine-kernel to the Airy kernel. 
In terms of the Airy function Ai{x), defined by 

A,{x) = — / rf^et^ +^"^, (42) 

which satisfies the differential equation 

A'^{x) = xAi{x), (43) 

one has 

A-(.,y) = ilMiM^iS(£)iiM. (44) 

x-y 

In (^^ we have used the scaling variables x and y proportional to N'^^^{X — Ac) . 

There are two ways to obtain the moments under consideration. The first one is 
to write as before 

I, =< |det(A„ - X)r >= / 1^^^^ n "<) (45) 

but in this case, there are three periodic structure due to three valleys of Airy func- 
tions, and the result is more complicated. It does not seem to be expressible as simple 
products of gamma-functions. However we can use a direct method starting with the 
expression 

Ik =< [det(Ac - X)f^ >= j^JzA\z)ei^""^^' (46) 

This representation is the edge limit A; — * of 



'2K 



rdu- A(z)A(u) 

ran, ix[z)ix[u) ^ 3 

^^ V 27rz Hi H/ - Ac + Az ^ ^ 
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The sums and products over / run from / = 1 to / = 2K. The dependence of F2K on 



is of order 



We may then use the standard orthogonal polynomial method. To the complex 
measure 



we associate the orthogonal polynomials pn defined as 

Pn{x) = + lower degree, 

and 

J dlipn{x)pmix) = hn5n,m 

The integral of is then simply 

/ = Klhohi ■ ■ ■ hx-i 



(4^ 



(49) 



(50) 



(51) 



Note that this looks similar to the partition function of a matrix model, but here it 
is the partition function of a K x K matrix model, instead of N x N {K is finite, 
since it is the order of the moment that we are considering, whereas goes to 
infinity). Therefore this is for any K a completely explicit expression of the moments 
at the edge. Those coefficients /i„ are expressible in terms of ratios of determinants 
constructed with the moments of the measure : 



hr. 



dr, 



d. 



(52) 



n-l 



with 



dn = det 



/ mo 

mi 



mi 
m2 



m 



n+l 



(53) 



V m„ m„+i • • • m2n ) 
in which the m„ are the moments of the measure. Those determinants are constants 
along anti-diagonal lines (Hankel determinants). Then the product hohi ■ ■ ■ hx-i is 



9 



reduced to a single determinant. For example, we have for K = 4 



hQhih2h^ = det 







-lC2 





^C^ \ 











2C2 









2C2 





V 




2C2 





-4Ci/ 



(54) 



with Ci = Ai{0) = 3-2/3/r(2/3), C2 = A^(0) = -3-^/71(1/3), since all the moments 
up to rriQ are easily expressible in terms of mo and mi alone. 
More generally we have 

m„ = J z^d^i = (-i)"(n -2){n-5){n-8)---An (55) 

where An = Ci for n=0 (modulo 3), and An = C2 for n=l (modulo 3) and An = 
for n=2 (modulo 3). The last parenthesis of the product in the r.h.s. of ( p3| ) is 
the rest of the division of n — 2 by 3. Then, dn is the determinant of a Hankel 
matrix, whose matrix elements in the first row are [< >,< z >, < >,■■'] = 
[Ci, -iC2, 0, iCi, 2C2, 0, -4Ci, 102^2, 0, -28iCi, -8OC2, 0, ■ ■ ■], and all the others are 
given by the Hankel rule. In this way we obtain sucessively, 

= Ci = 0.355028053 
hohi = Ci = 0.066987483 



2 



Y[hi = 2Cl + C'l = 0.010074161 


3 

Y[hi = -8(:7i(:73 - 3Cf = 0.001580882 


4 

l[hi = 12Cl + 2W\Cl = 0.000313095517 


5 

l[hi = -2160^6 - 1952C3C| - 432Cf = 0.000090756324 (56) 


Therefore for the edge problem we have found moments given by 's which are 
more complicated since = Ilo'^"^ hi. The result is explicit for any finite K, but we 
have not succeeded to continue it to non-integer K. The numerical values indicate a 
smooth curve in a logarithmic plot. 
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5 Finite N results 



We have derived in our previous paper the correlation functions of the character- 
istic polynomials under the form of a determinant. 



K 



i^i^(Ai,---,A/<) = (ndet(A;-X)) 



1 



A(Ai,---,Ax) 



det 



Pm(Ai) Pm+i(Ai) 
Pm{>^2) Pm+i(A2) 



Pm+k-i{^i) 



(57) 



Pm{^k) Pm+i{^k) ■■■ Pm+k-i{^i 
in which X is an M x M random matrix. 

The polynomial Pn{x) are the (monic) orthogonal polynomials, whose coefficients 
of highest degree are equal to unity 



Pn{.x) = + lowerdegree. 



(58) 



If we are concerned simply with the moments of the distribution of a single char- 
acteristic polynomial, we obtain from (|57|) 



/ix(A) = Fk{\ ■ ■ ■ , A) = ( [det(A - X)]^ ) 

Pm(A) Pm+i(A) 
(_l)i^(i^-i)/2 p',,{\) p'm+i(A) 

~ — 1<' ^ , . — det 

n£i'('!) 



Pm+k-iW 
p'm+k-iW 



Pm+k-iW 



(59) 



For the Gaussian distribution, 



with 



P(X) = ^exp-^TrX^ 



M = N -K, 



(60) 



(61) 



the polynomial Pn {x) are the Hermite polynomials {x) , defined with our normal- 
ization as 

1 V (jjJb 
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The integral representation 
HJx) -- 



AT" 



2iTT 



y{n+l) 



(63) 



over a contour which circles around the origin in the z-plane, turns out to be well 
adapted. 

Note that all these expressions are all valid for finite N. We may thereby recover 
readily several results that we have discussed in the previous sections. For instance, 
let us assume that M is an even number, and consider the center value A = (since 
the dependence in A is known to be contained entirely in the overall factor [p(A)]^^, 
as far as the coefficient •jk is concerned, it is sufficient to put simply A = 0). 

The Hermite polynomials Hn{x) vanish for odd n at x = 0. Similarly the odd 
derivatives of Hn{x) for even n, also vanish at x = 0. Hence, the elements of the 
determinant (^) are alternatively non-zero then zero. Thus the determinant is de- 
composed into a product of two determinants ; this is the exact phenomenon for N 
finite of the factorization of the symplectic and orthogonal determinants that we have 
seen earlier for large N. Since the matrix elements of (|59|) at A = are all expressed 
as derivatives of Hermite polynomials at the origin, it is possible to compute this 
determinant exactly for finite and arbitrary M and K . For the even M case. 



F2k{^) 



_l\K{2K-l) 



11«=3D0 



(/!) 



det 



H'iM 



H 



M+2 



(0) 



H 



{2K-2) 
M 



(0) 



{2K-2) 
M+2 



(0) 



X det 



M+3(0) 



-f^fl"_i_3(0) 



H 



M+1 (Uj iiM+3 (Uj 



(64) 



We denote each determinant as j(^)/A^^*^/2 and /(2)/jY-^Af/2 respectively, and 



2K-1 ,,• 
1=0 



(65) 
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The above two determinants are easily computed through the exphcit expressions for 
the if„(x)'s, 

ff.„(.) = i(-ir(2„ - 1).. f ••:'-"+'",:^' A(^r. (66) 



n 



m=0 V 2 ^ 2 



(i)(| + l)---(i + m- 1) mr 2 



1 / .^n.o^ , .^n^^ (-n)(-n + l)---(-n + m-l) 1 ^Nx\^ 



H2n+i(x) = (-l)"(2n + l)!!x y — ^ ^ , . 

(67) 

The two determinants contain overall products of factors of the form {2n — 1)!! ; once 
they are extracted one finds 

/(I) = C{M + 2K -3)\\{M + 2K ~5)\\---{M 
- r ^ f. T{M + 2l-2) 

/(2) = C(M + 2i^-l)!!(M + 2i^-3)!!---(M + l)!! 
^ 1 ^^ V{M + 2l) ^ 

= C ^(M+K-l) Hi -n/M , ,N J (69) 



with 



C = 2^^ n ^' (70) 
which is independent of M. In the large M limit, from the Stirling formula, we have 

,,s MK+K(K-1) MK K 

/{!) ~ CM 2^ -e'—2- (71) 

/ON „ MK+K(K+1) MK K 

/(^) ~ CM 5 e-—2- (72) 

It is remarkable that, even for finite M (M is the size of the random matrix), F2k{'^) 
for this Gaussian distribution, already exhibits the factor = H^o^ ^'/((-^ + 0' = 
[n^o^ ^']^/ nf5)"^ which is known to be universal in the large M-limit. It is indeed 
obtained from the product of the factor C and 1/(11^^"^ ^0 in (|65D . This means that 
at each order of the expansion, we keep this universal factor for -F2i<'(0). In the 
large N limit (M = N - K), F2i<-(0) becomes 

F,^(0) ^ (2iV)''^e-"^^ n'^;^^- (73) 
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In the previous paper, we have derived F2k{^), in the large N hmit, as 

FMK ■ ■ ■ , A) ^ (27rp(A)iV)^^e-^^ Jl (74) 

At the band center, A = 0, the density of state is p(0) = ^ for the Gaussian distribu- 
tion. Therefore, ( [73D is indeed consistent with (^). 

It may be interesting to note that one of the factors of (0), namely llf=o~\l^-), 
appears in F2k{^) in (^) . This factor, a product of gamma-functions, remains for 
any set of orthogonal polynomials, since it stands in the front of the determinant of 



The factor e~^^ is cancelled by the normalization ||^. For A 7^ 0, we have 
evaluated F2k{,^)- We have here considered the finite N case to see the universal 
factor 77^-. 

One can recover again the Airy limit by the use of eg. (|59D . We use once more the 
properties of the Hermite polynomials such as 

H'^{x) = nH{x) (75) 

and their explicit integral representation 

H^{x) = ^iV^e^"' r dss^'e-^''-'^'"'. (76) 



00 



/27r 

We set n = 5 + N, and after exponentiation, we have 



Hn{x) = ^A^ieT-' r ds/e-^^(^) (77) 



^277 



00 



where f{s) = -|- isx — log s. The saddle points are degenerate at the edge x = 2. 
The vicinity of this point is blown out through a change of variables, with a scaling 
ansatz, 

x = 2 + N-^'y (78) 

and 

s = -i + N-f^z (79) 
14 



If one expands /(s) up to order 2;^, one sees that in the proper scahng choice a = 2/3 
and f3 = 1/3, one recovers the Airy hmit which governs the properties of the system 
in a neighbourhood of size gner's semi-circle. Then, the 

integral becomes 



This is indeed the Airy function Ai{z) of ( ^2]) 



Hn+s{x) = V2^e^^{-iYAi{{x - 2)A^3) 



We now consider all the Aj = 2, and the determinant (|59D becomes in the large N 
limit a determinant of Airy functions. If we replace Hm+2K-i at the right-up corner 
of the determinant by the Airy function Ai{0), the other matrix elements become 
derivatives of the Airy function, since there is a the recursion relation (^). For 
example, in the K = 1 case, we have 



det 



^m(2) ^m+i(2) 



det 



Am A{0) 



Then, we find in the large N limit, with 



F2k{2) 



nt 



2K-1 



=0 



det 



2 

M+1 



= M-K, 



m) MO) 
^r(o) Am 



12) 



^3) 



The above determinant was discussed earlier. Note the factor l/nz=o ^' front. 



6 Derivative moments 



The same techniques may also be used if one is interested in the moments of the 
D-th derivatives (D = 1,2,...) of the characteristic polynomials. Let us consider for 
instance 



<D). 



2K qD 

2Ki^u---,^2K)=< n 7^det(Ai-X)> 

l=3Dl '^'^i 



15 



From (^9]), one sees immediately that it has also the form of a determinant : 



A(Ai,---,A2i^) 



det 



pS^(Ai) pi?|i(Ai) 
Pm(A2)(^) A?|i(A2) 



(Ai: 

PM+2i^-l(A2; 



When all the Aj's are equal, we have 



4?(A,---,A) = {[^dei{X-X)f^) 



det 



p1?(A) A?|i(A) 



(D+2K-1), 

Pm ' 



(A) Pi?:r-'^(A) ■■■ p'mI^kZ^'W 



Pm+2K-1 (A) 

(A) 



^6 



If we set A = it may be again decomposed into a product of two determinants . 
Let us assume for definiteness that both M and D are even. Then, we have 



J(^) = det 



H 



(D+2A'-2) 
M 



(0) H 



(D+2K-2) 
M+2 



(0) 



(87) 



= det 



Using the explicit expressions for the Hermite polynomials , we can compute these 
determinants. We find for arbitrary M,D and K, 

1 o^ 1 



]S[K(M-D) 



9) 



/(I) 



{M + 2K- 3)!!(M + 2K- 5)!! ■ ■ ■ (M - 1)!! 
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X n[(Y + 0(Y + ^-i)---(Y-^ + ^ + i)] 



DK+K(K-1) ^ ^ 



X 



2 -2 n ^- (90) 



1=0 



/(2) = (M + 2K- l)!!(M + 2K-3)!!---(M+ 1)!! 

X n[(Y + ^ + i)(Y + 0---(Y-^ + ^ + 2)] 

X 2 5^ — - n ^! (91) 

z=o 

(One may easily check these results for D = M, since the matrix elements below 
the diagonal vanish, i.e. the determinants are then simply given by the product of 
the diagonal elements, Ilj^oiM + 2l)\ which agrees with (pO|) . When D = 0, it reduces 
to the previous expression (H). is obtained from J(^) by the shift M ^ M + 2.) 

In the large N limit, we have 

f,',?'(0) i2Nf'^-""'e-'"'-^ n (92) 

Hence, for this derivative moments at finite M, again the universal factor 7^ is 
present, and it persists of course in the large N limit. 

These results lead to the conjecture that the average of the moment of derivatives 
of the Riemann zeta-function along the critical line 

also have this universal factor 'Jk- 

7 External source 

We now consider the case in which the external source matrix A is coupled to the 
random matrix X. The measure of the random matrix X is 

rf^(X) = 1 -^T-^^'+^^rf^'x (94) 

Zj 
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The eigenvalues of the matrix A are denoted by ai, i = 1, N. In such cases, 
the standard orthogonal polynomial method cannot be used. However, the n-point 
correlation functions p(Ai,---,A„) have been found to be described again by the 
determinant of a kernel ; from there the level spacing probability p{s) has been also 
investigated |0. If we specialize to a source which has two opposite eigenvalues, 
namely = +a for i = 1, ■ ■ ■ , N/2 and = —a for i = N/2 + 1, ■ ■ ■ , N, one finds 
a support for the eigenvalues made of two disconnected segments for a > 1. If one 
tunes the external source so that a = 1, i.e. = ±1, the gap between the two 
segments closes and the spectrum consists of a single segment for a < 1. We want 
to investigate here the critical point a = 1 which gives rise to yet another class of 
universality. 

The moments F2k{^, ■ ■ ■ , A) at A = at this closing gap point may turn out to 
have interesting applications. 

Since X and A are Hermitian matrices, we write 

TrXA = TrU^^XoUAo (95) 

where Xq = diag{xi, ■ ■ ■ , xm), Aq = diag{ai, ■ ■ ■ , Qm), and U belongs to the unitary 
group. The integration over this unitary group U is well known from the work 



of Harish-Chandra, Itzykson-Zuber [p!8| , [T^, and this is the starting point of the 
formulae found in For instance the n-point correlation functions are given by the 
determinant of the n x n matrices made with the kernel K{Xi, Xj) with 

J-oD 271 J 2ni fJi u — ai u — it 
where the contour encloses all the a^'s. 

However we may proceed without that here and compute the correlation functions 
of the characteristic polynomials directly. Indeed 

K 



Fk{\i,---Ak) = <ndet(A„-X)> 

= ^ [dxfl det(A, - x)e-TTtx'+NT.xA (g^^ 
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a=l 



In the above equation, the random matrix X is assumed to be an M x M matrix, 
with M = — i^, as before. When K = 1, this gives a polynomial, which was 



investigated before [17 



The explicit integration over the unitary group ||T8|, |1^], leads to 

J fJi A(a)A(A) 
where A(xi, ■ ■ ■ , xm] Ai, ■ ■ ■ , Xk) is the Van der Monde determinant (M + K) x (M + 
K) made with the x's and the A's. This determinant may be replaced by a deter- 
minant of (monic) polynomials , and we choose the Hermite polynomials defined in 
It is then straightforward to verify that 

r i7„(x)e-T-'+^-dx = a^eT'^V^ (99) 

J — oo V ' 

Therefore we can explicitely integrate over the M variables Xj's in (^8|) and one obtains 

= A(a)A(A) 

1 ••• 1 ^o(Ai) ... HoiXK) 

xdet : ■•. : : ••. ; . (100) 

Let us first check that in the limit of a vanishing source in which all the Oj — 0, we 
do recover the previous formula (|57|). The column which depends upon is expanded 
in Taylor series around Oi, and subtracting the successive columns, we obtain, after 
factoring the Van der Monde determinant A (a) which cancels the denominator, a 
vanishing upper triangle (up to the M-th column) , ones on the diagonal and powers 
of the tti's below the diagonal. We can now let the a^'s go to zero and we are left 
with the K X K of (pTf). (In we gave a different derivation of this same formula). 

If we return to an arbitrary non vanishing external source, we may proceed by 
returning to ( [L00| ) and define Gxipi, ■ ■ ■ , 6/^), 

= ^efi:^? (101) 
A(a) 
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We may now recover Fx by taking the Fourier transform of Gxiibi, ■ ■ ■ , ibx] 

^ dbi , 1 

i=\ 

Therefore, we obtain the following explicit formula, 



/ GK{^h,, ■ ■ ■ , ^hK)e-'^T.^^'^ ft ^ = (^)''A(A)i^i.(Ai, ■ ■ ■ , A^)e-f (102) 



M K 

X 



1=1 j=i 

(103) 



Note that we could replace in the integrand of ( |103| ) det(e ^^^^^i') by the diagonal 



term e~* ' ' and cancel the i^'! in the denominator. Again we can examine the 
limit of this formula when the external source goes to zero, and putting all Aj = A, 
we obtain 

(we have considered F2K instead of Fx in order to compare with our previous results). 
In the large N limit, we exponentiate bf'^ , (M = N — K), and look for the saddle 
points which are the roots of the equation 6^ + iXb — 1 = ; let us call the two 
roots b'^ and b^. The difference |6+ — b^\= 27rp(A). The leading saddle-point for the 
bi's,{l = 1, ■ ■ ■ , 2K), is obtained by choosing half of them equal to 6^, and b^ for the 
another half. The following Gaussian integral with a Van der Monde determinant, 

1=1 i<j J (iV J j 2 

where /" is the second derivative of / at the saddle-point, allows us to complete the 
calculation. Integrating then around the saddle-points 6"*" and 6~, and keeping in 
mind the combinatorial factor ^f^| , which is the number of choices of K b'^ and K 
b~ among the 2K bi^s, we recover precisely our previous result, 

e-^^^WF2i,(A) = (27rArp(A))^'e-^^7^ (106) 
20 



where 7^ = {Uto' ^!)V(n?fo'' ^0 = (H/Io' ^0/ Uto\K + /)!, and V{X) = f . 

From the expression ( |103| ), it is also easy to obtain the moments at the critical 
point corresponding to the closure of the gap : 

FK{0) = ^e— J Jl^{l + bf)-A\b)e--^i-r'' (107) 

Note that this expression is exact for finite N. In the large N limit, we exponentiate 
the logarithmic term and expand the exponent about bi up to order bf term. The 
critical point, is precisely the point at which the coeffeicient of the quadratic term bf 
vanishes. We then have 

„ . , A^"'^ JVM f 1 r db] N f,4 . o X 

Fk{0) = -j^e^ jYl-le-TE^^.^A'ib) (108) 

As in all the cases which appeared in the previous sections, this integral is expressed 
by a Hankel determinant, in which the matrix elements are r(^^j^). The determinant 
is 



/ = det 



r(i) r(f) 
r(|) r(f; 
r(|) r(f) 



(109) 



Note that we have considered the case a; = ±1 case, but the formulae are explicit 
for any spectrum of the source and they could be easily used to study for instance 
multi-critical situations which were discussed in 0. 
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